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Abstract 

We consider ideal chains in a hypercubic lattice d > 3, with a fixed ratio m of self-intersection 
per monomer. Despite the simplicity of the geometrical constraint, this model shows some inter- 
esting properties, such as a collapse transition for a critical value m c . Numerical simulations show 
a Self- Avoiding- Walk-like behavior for m < m c , and a compact cluster configuration for m > m c . 
The collapse seems to show the same characteristics as the canonical thermodynamical models for 
the coil-globule transition. 
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I. INTRODUCTION 



Polymer models have been the subject of extensive theoretical and numerical studies. 
Most of these models are based on ideal flexible chains, with the addition of various kinds of 
interactions between monomers, to include some non-trivial properties of real systems [HE]- 

Very important is the excluded volume effect, which significantly modifies the fractal 
properties (the chains behave like Self- Avoiding Walks, SAW); and the coil-globule (CG) 
transition, in which a flexible chain collapses from an extended coil to a liquid-like globule 



Most of the thermodynamic models that show the CG transition usually consider a com- 
petition between interactions of different geometrical nature. For instance, one can consider 
a lattice random chain with repulsive on-site interactions and attractive nearest-neighbors 
links. The transition arises from the competition of these interactions. 

Let uj = {Sq, Si, ... ,S n } be an n— step Simple Random Walk (SRW) on the Z d lattice 
(Si are lattice vectors), and define the number of visits to each site x as ip x = S (x — Si). 
A canonical model, incorporating excluded volume and a CG transition, is described by the 
Hamiltonian 



where (x, y) indicates nearest-neighbors (taking e± = will lead to the Domb- Joyce model, 
see |1]). 

A few years ago an interesting thermodynamic model that shows a transition with only 
on-site interactions was proposed [5]. However, a competition between two geometrically 
different constraints is still present since a self-avoidance is incorporated through restricting 
the maximal number of visits per site. 

To our knowledge all the thermodynamic models proposed so far, having only one kind 
of short-range interaction, do not allow for the CG transition in the Boltzmann parameter. 

In this work we will show how excluded volume effects, the CG transition and liquid-like 
clusters can be obtained by imposing a single global geometric constraint. We consider 
ideal chains in a hypercubic lattice Z d , d > 3, with a fixed ratio m of self-intersections per 
monomer. 

This model shows a CG transition for a critical ratio m c : numerical simulations (obtained 
by a standard implementation of the Pruned-Enriched Rosenbluth Method, PERM, see j6]- 



[2112]. 
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9j) indicate a SAW-like behavior for m < m c (for a review about SAW see [T0HT2] ) and a 
cluster configuration for m > m c . In addition the model is amongst the simplest with a 
crossover from SAW to cluster behaviour. Our focus will be on d > 3 lattices because in 
these lattices, as we shall see, the transition occurs at non-trivial values of m c . 

II. MODEL DEFINITION 

Consider an ideal chain u of n steps on Z d (without loss of generality we take So = 0). 
We call R[u] the range of u (number of distinct lattice sites visited by the path), and 
M [uj] = (n + 1) — R[w] the number of self-intersections. 

Let P n (M) be the fraction of SRW, of length n, with exactly M £ [0, n - 1] self- 
intersections. We introduce mean value (M n ) and variance (AM 2 ) = (M%) — (M n ) 2 of 
the distribution P n (M). It is well known (see [4j) that for SRW in d > 3, d ^ 4, 



where d is the probability that an infinite length walk contains its starting site at least 
twice (for numerical values of d, see [13]), and A^ is exactly known (for d — 4 the main 
fluctuation is actually A 4 log(n)) jH [13]. Concerning the variance, Jain and Pruitt have 
shown that (AM^) cx nlog (n) for d = 3, and (AM^) oc n for d > 4 |U [H]. They have also 
shown that for d > 2 



is normally distributed, from which it follows that P n (M) is peaked around its mean value 
for long walks. 

We are interested in chains in which the rate of intersections per monomer is fixed at a 
certain value m £ [0, 1] in the thermodynamic limit. We define the ensemble fl m of n— step 
walks with exactly M — [m (n — 1)J intersections, and Q, the ensemble of all n— step SRW: 
we call (*) m the average on Q m , and (■) that of Q. 

For large n we can approximate \m (n — 1)J ~ mn. Henceforth we will work under this 
approximation. We introduce the fraction of m— intersection rate walks as P n (m). From 
the properties of P n (M) it follows that P n (m) is peaked around (m) = Cd when n — > oo. 

Let (S%)m De the mean square end-to-end distance on Q m : from our simulations we find 
that the relation (S^) m oc n 2ud(jri ^ holds for any dimension considered, with the exponent 



(M n ) =C d n-A d n 2 -? + O {!) , 



(2) 



£= hm (M-(M n ))/(AM 2 ) 



(3) 
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Vd (m) dependent on m and d. Again, for d > 3, the simulations show the existence of a 
critical value m c £ (0, 1) beyond which the chains collapse into a compact liquid-like globule, 
with Vd (m) = l/d. If instead m < m c we have the SAW-like behavior (Fig. (fil) and Q). As 
preliminary observation we can state that, since for m = we have the Self- Avoiding Walk, 
then (S%)o — DdT? Vd [yd is the correlation length exponent for the SAW, see jTD]) ■ On the 
other hand, from Eq. ^ it follows that the relation 

f " de (S 2 n ) Cd+ ePn {C d + e) ~ n, (4) 
should hold in the large n limit and b n = o (1). This is clearly confirmed by our simulations. 



III. GENERAL RESULTS 



We will briefly discuss the recurrent cases d = 1,2. For d = 1 we have a rather simple 
situation: in a linear lattice chain the end-to-end distance is proportional to the range, hence 
proportional to n, Vm < 1. A unidimensional SRW is recurrent [C\ = 1, see [I]), therefore 
no drop of V\ (m) is expected for m < 1. At d — 2 the picture is conceptually similar. 
We find v<i (m) = (data in Fig. Q), and since the square lattice SRW is still recurrent 
(C*2 = 1, from |3]) from Eq. Q we expect that this relation holds for any m < 1 (see later 
discussion about the case d — 3). 

The most important case is clearly c? = 3: from simulations we find that z/ 3 (m < m c ) = 
u 3 , and z/ 3 (m > m c ) = 1/3. The last statement needs some attention, we point out that 
U3 (in > m c ) = 1/3 is reached only after a long crossover. Fig. shows v vs rT x M and 
C3M/ (M n ): the former does not consider the finite size of the system, while the latter takes 
into account fluctuations of the walk mean support (range) (R n ) at finite n (see |15|). This 
enhancement allows a significant improvement of the accuracy for d = 3 (while is useless in 
higher dimensions, since, by Eq. (|2]), the range converges rapidly to its asymptotic behavior). 

From Eq. ^ we can locate m c for d = 3. If (S%} m nas the SAW-like behavior (with 
z/3 > 1/2, see [TT]) for m < m c , and cluster-like for m > m c , then the drop of z/3 (m) must 
lie within the range of integration of Eq. (|4]). It follows that m c = C3 for n — > 00. As for 
the exponent at the critical point, Eq. ^ tells us that Vd (m — > Cd) = 1/2, from which it 
follows that z/ 3 (m c ) = 1/2: this observation is supported by Fig. Q and ([2]), where we see 
that 2z/ 3 (m) passes through 1 at the expected critical point m c = C3. 
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Tabella I: Summary of predictions for Vd (m): for d = 5, 6 we where unable to look at the critical 
behavior, while the critical point m c is evaluated by the drop zone of fd(m). The I.e. indicates 
SAW logarithmic correction in d = 4. 

This latter fact is of some importance, since the transition for d = 3 would show the 
same behavior as that described by the Hamiltonian in Eq. Q. If the two models belong to 
the same universality class this would be very interesting, since in our model the transition 
arises from the necessity of maximizing the configurational entropy, without the need for any 
further interactions. Indeed, at the critical point m c there is a radical change in the optimal 
strategy to achieve the global constraint m (which is actually a long-range correlation): for 
m < m c , the best way to change the ratio m is to compress (or expand) the chain locally, 
keeping the SAW-like fractal structure. In this situation the monomers intersect (on average) 
only within a certain distance along the chain. For m > m c instead, it becomes entropically 
convenient to assume a compact configuration, also allowing intersections between monomers 
very far apart in terms of position along the chain. 

We also studied higher dimensions, the results are consistent with a CG transition of 
the same kind. For d = 4 we have 1^4 (m > m c ) = 1/4 and z/4 (m < m c ) = 1/2, but the 
random-coil behavior shows logarithmic corrections. Simulations fit the conjecture that 
(5^)m 00 ^log (n) 1 ^, as in the SAW case [10] . If this is true, from Eq. we should again 
find m c = C4. 

As expected, for d > 5, m < m c , we find the mean field behavior 1/4 (m) = 1/2 (as for the 
SRW), and Vd (m) = 1/d for m > m c . Since for d > 5 Ud (m < m c ) is the same as for the 
SRW, we neither use Eq. Q to locate m c nor find Ud(m c ), but our simulations (extremes 
of the drop zone of i/a (w), d — 5, 6 in Fig. Q) strongly suggest that m c > Cd- All results 
and conjectures about the behavior of {S%)m have been summarised in Table I. 

We would like to point out that all results for d > 3, m > m c , are consistent with a 
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remarkable exact work by M. van den Berg, E. Bolthausen and F. den Hollander Ql6j) on 
the moderate deviations for the volume of a Wiener Sausage (WS), which is a neighborhood 
of the trace of a standard Brownian motion up to a time t, given by taking all points within 
a fixed distance a of Brownian motion (essentially a continuous version of our model, see 
|17j). Let 77 (t), t > be the standard Brownian motion in M. d starting at the origin. The 
WS W a (t) with radius a is the process defined by 

W a (t)= |J B a ( V (s)), (5) 

0<s<t 

where B a (x) is the open ball with radius a around x G lR d . This paper considers the 
probability of having a Wiener Sausage W a (t) of volume \W a (t)\ < bt, b G [0, K a ] (« a t is 
the long time behavior of (|W"°(t)|)), showing that there exists a critical value b c below 
which the sausage is supposed to collapse in a swiss-cheese like compact configuration (a 
non-percolating cluster with random holes of size (9(1)). They rigorously showed that 
b c G (0, Ka) for d > 5 only, while for d = 3, 4 the transition is at b c = K a exactly. 



IV. ORDER PARAMETER 

Typically the order parameter considered for a CG transition is n~ 1 (S'^) / 3 (with (-)p being 
the thermodynamic average); we propose here a slightly different parameter: 



Qd 



(m) = lim y/(S*)J(S*) . (6) 



This function vanishes beyond the critical point, similarly to the magnetization for spin 
systems. If the analogy with the CG transition from Eq. (jlj holds, we should expect a 
second order transition, at least for d — 3. 

Our simulations (Fig. Q) do not allow us to clarify if the transition is continuous at 
m c . However, for d = 3 we can get some insights by a mean field analysis of the Stanley 
Model (SM, see [18J), a model of correlated random walks defined by the following partition 
function: 

n 

Z n (P) oc^Pn (M) e-? M . (7) 



M=0 



Given that variable £ in Eq. ^ is normally-distributed, it's easy to show that for small (3 

lim n-\M n ) p ~C 3 - 5m (/?) , (8) 



where (M n )p is the mean number of intersections in the SM at temperature (3, and Sm (/?) 
vanishes for /3 — >■ 0. From Eq. ^ we have a Gaussian shape P n (M ) near M = C^n, with 
(AM^) oc nlog (n). By replacing M — > mn in Eq. ^ we find 

Z n (/?) oc f fee-f^fAMS- 1 ^] (9) 

for small /3, with Sm = (C3 — m) ~ 0. Whereas there is a positive constant a such that 

(AM%) > a^n for large n, performing a saddle point integration on Eq. ^ with n — > 00 we 

obtain 5m > Co/3 (for small /3 and cq > 0). 

From simple mean field arguments (see [21]) it is reasonable to assume that the SM can 

be approximately described by the Flory theory. Consider the following Flory energy 

2 

Tp (r) ~ (3n 2 r- d + T — - [d - 1) log (r) + A (n) , (10) 

where r = |i? n | is the end-to-end distance, and A (n) is a function independent of r [in]. We 
minimize the functional under the assumption r ~ (3 e n u , thus 

rf/3 i-0(i+<V-Ki+d) + (d - 1) (i' e rT v = f3 e n u -\ (11) 

When ud<2we find z/ = 3/ (d + 2) and 9 = 1/ (d + 2), otherwise v = 1/2 and 6 = 0. When 
d = 3 this theory predicts 8 = 1/5. A different estimation method uses path integrals, results 
are in agreement [22]. 

Flory theory for SM where d = 3 predicts (S%)p oc /3 2//5 n 6//5 for small (3: invert- 
ing Sm(f3), and substituting in the expression for (5 r 2 t )/3, we find £3 (m) < boSm 1 ^ 5 , with 
Sm = (C 3 - m) ~ 0, 6 > 0. 

In the case d = 4 Flory theory cannot be used, since the logarithmic correction to the 
mean field behavior are not captured. Note that the d > 5 theory predicts 9 = 0, in 
agreement with the prediction m c > Cd for d > 5. 

The exponents found are certainly incorrect, since z/3 7^ 3/5, but a power law behavior 
£> 3 (m) oc Sm e (with a likely logarithmic correction) is reasonable, suggesting that in d = 3 
the transition is of the second order, with m c = C3. 

V. CONCLUSION AND OUTLOOK 

In this paper we introduced a new athermal model of interacting random walks, which 
shows a CG transition for a critical ratio between the range and the number of monomers. 
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The transition seems to show the same characteristics as that seen in canonical models. 
However, the relationship between our microcanonical model and those canonical for the CG 
transition is non-trivial (for a general discussion on the differences between microcanonical 
and canonical ensembles see |23|). 

The Hamiltonian in Eq. ([!]), as well as the interaction described in reference |5], operate 
a selection on the ensemble Q whose mechanisms are not easily connected to the ensemble 
Q m . Although this is a key issue, we do not discuss it here: a work dedicated to this 
subject is currently under preparation, where we will present a detailed study of the P n (m) 
distribution and its relation to some thermodynamic models. 

Apart from more focused implementation of the simulations here presented, other issues 
of interest could affect the connectivity properties of collapsed clusters (an example is found 
in [16], collapsed WS clusters should be non-percolating, with O (n) holes sized O (1) dis- 
tributed inside the range: is reasonable to expect the same behavior for our model when 
m > m c ). 

In general, this model certainly deserves attention since, in our opinion, it could lead to 
substantial improvements in understanding the geometry of CG transitions, as well as the 
crossover between the SAW and the SRW. 
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Figure 1: Exponents Vd{m) vs m/Cd from simulations at various lattice dimensions. The graph 
shows the behavior of (m), obtained from fits of log ((S^)m) vs log (n). The range for all fits is 
n £ [0.8 • iim, um\- Maximal length for simulated walks was: um = 2 • 10 3 for d = 2, um = 10 3 for 
d = 3, um = 10 3 for d = 4, um = 0.5 • 10 3 for d = 5, and um = 0.3 • 10 3 for d = 6. Moreover, for 
d = 3, the log ((5^) m ) has been achieved by using m = C^Mj (M n ), with an n dependence, to take 
into account the finite size of the chains (see Fig. (|2])). For d > 3, the finite size effects are much 
weaker and this correction is not necessary. The picture for d = 4, m < m c shows an exponent 
slightly larger than 1 (a correction of order O (lO -2 )), consistent with the expected logarithmic 
correction (for chains of length 10 3 an exponent of 2v ~ 1.070 would be predicted). In the n — > oo 
limit, Vd (m) is expected to be a step function (see Fig. ([3])). The present graph shows the step at 
m c : for d = 3, 4 we expect m c = Cd, for d > 5 simulations suggest m c > Cd.. 
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Figure 2: Exponent (m) vs m, a comparison of simulations with and without finite size correc- 
tions. The chart, from fits in the form log ((S^)m.) vs log (n), compares m = M/n (empty circles) 
with m = C%M / (M n ) (full circles), (M n ) is obtained from Eq. §2^. Essentially, M is rescaled with 
the average value of the support (range) {R n ) = n — (M n ): this procedure leads to a significant 
improvement of the accuracy, at least for m > C3 (vertical dotted line). Given that the two curves 
in the limit n — > 00 must converge, it is clear that for hm = 10 3 we are still far from the asymptotic 
regime. This fact is emphasized in Fig. ([3j in which we can see how the drop band for d = 3 
vanishes very slowly compared to d > 3. 
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Figure 3: Slope of Vd (m) in the drop band. The graph shows (in log-log scale) the maximum slope 
Gd (n) = — [dm^d ( m )] m i n °f u d ( m ) f° r different lengths n, which is an estimator of the drop band 
width (proportional to G^ 1 (n)) in which Vdijn) falls from v& to 1/d. Gd (n) seems to increase with 
a power law, supporting the idea that Vd (m) is a step function in the limit n — > oo. Dotted lines 
are power law B ■ n a fits to the data: a = 0.29 ± 0.01 for d = 3 (without the finite size correction 
C 3 /(M n )), a = 0.67±0.03 for d = 4, a = 0.83±0.05 for d = 5, and a = 1.02±0.07 for d = 6. From 
this data we can notice that the width of the drop band decreases faster in n as lattice dimension 
d increases. 
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Figure 4: Order parameter Qd( m ) vs m /Cd from simulations at various lattice dimensions. The 
graph shows (S T 2 t ) m /(^) o: „ = 10 3 for d = 3, n = 2 • 10 3 for d = 4, n = 10 3 for d = 5 and 
n = 0.5 • 10 3 for d = 6. The vertical dotted line is the expected transition point for d = 3, 4 as 
n —7- oo. The asymptotic behavior near m c is only hinted, since our PERM implementation, which 
requires us to store the microcanonical density of states at many n, has not allowed us to simulate 
longer chains. Better pictures could surely be achieved, by more focused techniques, for every 
dimensions except d = 4. For d = 4 and m < C4 we expect n~ 1 {S^) m to be only C(log (n) 1//4 ), 
implying the asymptotic behavior would be difficult to demonstrate clearly even by simulation for 
very large chains. 
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